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I. INTRODUCTION 



Among several methods in nuclear theory, the density functional method provides one of the most widely applicable 
treatments of many-nucleon systems. For example, the nucleus is a many-nuclcon system, where more than 300 stable 
nuclei are already known. All atoms contain nuclei, and the chemical properties are identified by the number of 
protons included in the nucleus (cf. chemical elements: H, He, Li, • • ■ ). It follows that the nucleus is one of the most 
important ingredients of our universe. 

The theory based on the nucleonic degrees of freedom is presented, where the nucleon is not an elementary particle 
but consists of quarks and gluons. The validity of this theory is based on the fact that nucleons are quite stable 
quantum entities that can be regarded as effective physical units. Since nucleons are fermions, there are several 
common features relating this many-fermion system with other physical systems such as many-electron systems and 
so on. However, there are several specific features in nuclear energy density functional. Indeed, nucleons, which have 
isospin and spin degrees of freedom, interact by two completely different forces: nuclear and Coulomb forces. In 
addition it is worth noting that reactions between nuclei are not similar to reactions of the other physical systems 
to a large extent. Depending on the isospin degrees of freedom, there are two kinds of nucleons; i.e., protons and 
neutrons. While we understand the Coulomb force well, much is not known about the nuclear force. Indeed, although 
the quantum chromodynamics (QCD) Lagrangian has already been established, its connection to the nuclear force is 
still developing. It makes many-nucleon system research quite difficult and also fascinating. 

From a scientific point of view, many-nucleon research is associated with clarifying the origin, existence, structure 
and reaction of chemical elements, where the origin of elements heavier than iron has not been understood well. It 
is an attempt to understand the time evolution of our universe with respect to the constituent chemical elements; 
"How and where were all the chemical elements created and why do they exist as they are ?" . Chemical elements 
which do not naturally exist on earth are called superheavy elements. These elements are artificially synthesized in 
the laboratory. Such a superheavy element research is related with clarifying the existence limit of chemical elements. 

Finite-body quantum systems are often investigated in the research of many-nucleon systems. Although the spherical 
shape is expected to be energetically favoured if the system is governed simply by surface tension, the nucleus has 
experimentally been shown to have several shapes: e.g., spherical, prolate, oblate shapes and so on. Such a research 
is associated with the structure of nuclei. 

It is reasonable to have a unified theoretical framework describing both stationary and non-stationary states. The 
nuclear density functional theory is a possible candidate. The basic equation of the many-nucleon system is the 
Schrodinger equation containing the Hamiltonian. Based on the independent particle motion, a possible form of the 
Hamiltonian is generally provided by 



The first term of the right hand side arises from the kinetic energy, and the other terms in the right hand side from 
the interaction energy. For the interaction part, the density functional theory in many-nucleon system is a theory 
describing all the interaction by several densities (for many-electron systems, see Refs. As is already mentioned, 

the interaction part of the Hamiltonian is not perfectly known as far as the many-nuclcon system is concerned, so that 
the one important task in nuclear density functional theory is to find out the ultimate interaction, which can describe 
all the phenomena in the many-nucleon system. This point is indeed different from studying many-electron systems. 
The process of finding the ultimate effective Hamiltonian mainly consists of two steps; first, to find an appropriate 
functional form of the effective Hamiltonian; second, to find its best parameter sets. 

In this chapter much attention is paid to show a method of deriving the effective interaction in many-nucleon 
systems. What is presented in this chapter is a kind of modelling: the modelling of interacting femtometer-scale 
fermions. The nucleon-nucleon interaction is decomposed in Sec. [TTJ The effective Hamiltonian arising from the 
nuclear force is discussed in Sec. Mil starting from the zero-range nucleon-nucleon interaction (Sec. IIII A[) . we show 
a procedure of obtaining the Hamiltonian density (Sec. IIII B[) : the effective Hamiltonian is obtained by applying 
the variational principle fScc lIII Cj) . The effective Hamiltonian arising from additional forces is briefly discussed in 
Sec. IIVI In particular we provide a mathematically rigorous treatment of applying the variational principle, where 
infinite-dimensional Hilbcrt/Banach spaces are considered. For the mathematics used in this chapter, refer to the 
textbooks of functional analysis such as Ref. [lij . 
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II. NUCLEON-NUCLEON INTERACTION 



There are two kinds of nucleons (isospin degree of freedom r): protons and neutrons, which are fermions. In 
addition the nucleon has the spin degree of freedom. Therefore we arc interested in the many-fermion system with 
spin and isospin degrees of freedom. 

Nucleons interact by two kinds of forces: nucleon-nucleon interaction V is represented by the sum of the nuclear 
potential Vk, the Coulomb potential Vq and the pairing potential V pa i r : 

V = V K + V c + V pair . (1) 

Only protons interact by the Coulomb force. The two kinds of nucleons form a bound system called atomic nucleus. 
The atomic nucleus is a finite quantum system containing Z protons and N neutrons, where A = Z + N is called 
mass number. One of the goals of theory is to determine for which combinations Z, N can exist. 



III. NUCLEAR FORCE 



Let us begin with the primitive picture of the nuclear force using the potential description. The potential of nuclear 
force depends on the positions, momenta, spins (cr, = ±i) and isospins (t-£ = ±i) of the two nucleons: 

Vk = VK{Ti,Tj,Pi,Pj,<Tu(Tj,Ti,Tj), 

where i and j are the indices identifying the two nucleons. It is rational to determine interactions obeying invariance, 
because the physical law must be invariant regardless of observers. The following invariance are required for Vk- 
translational invariance, Galilean invariance, rotational invariance, isospin invariance, parity invariance, and time- 
reversal invariance (for each invariance, see Sec. 7.1.1 of Ref. |H). It is worth noting here that the radial dependence 
of the function Vk cannot be deduced from invariance principles. Among several attempts of determining the radial 
dependence, H. Yukawa proposed the Yukawa potential: 

Vy(r) = 

based on meson field theory, where 1 j \i is the Compton wavelength of the pion. 

Historically the components of the nuclear force are represented using the identity operator, (<Ji ■ (Tj), (t^ • Tj) and 
((Ti ■ <Tj)(Ti ■ Tj) multiplied by the spin- and isospin- independent ingredient. In the traditional formulation, such a 
nucleon-nucleon interaction is represented using several exchange operators: 

V K = V w (r) + V M {r)P r + V B (r)P a - V H (r)P T 
= Vw(r) + Vu(r)P r + V B {r)P a + Vn(r)P r P a , 

where P r , P a = i(l + cno-j) and P T = |(1 + TjTj) are exchange operators of coordinates, spins and isospins. In 
particular the fermionic relation P r P a P T = — 1 is used to replace the last term, and everyterm is represented without 
P T (for the property of fermions, refer to textbooks of quantum physics such as Ref. Q)- The indices of the first, 
second, third and fourth terms stand for Wigner, Majorana, Bartlett and Hciscnbcrg, respectively. Another important 
ingredient is the tensor force represented by (r^ • (Ti)(ri ■ (Tj), where the tensor operator is 

3{rj ■ aAjn ■ (Tj) 
bi * ~ uTja ' * ' i)~ 

The aim of this chapter is to have a density functional representation for the nuclear force (for a whole process, see 
Fig. 1), which includes the exchange properties represented by (<Xj • (Tj), (t^ • Tj), {(Ti ■ o"j)(tj ■ Tj) and (r,; ■ <Xj)(rj • (Tj). 
The interaction described by spin and isospin degrees of freedom is specific to nuclear physics. 
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FIG. 1: A process of obtaining the effective interaction. Ev K denotes the interaction energy, Vk is the zero-range interaction, 
denotes a trial function, and Hy K (r) is the interaction part of the Hamiltonian density. V e g, which is obtained by calculating the 
Gateaux differential, is the effective interaction in the density functional formalism. The zero-range interaction Vk is presented 
in Sec. IIII Al the energy and the Hamiltonian density are calculated in Sec. IIII Bl and the effective interaction is obtained in 
Sec lHTCl 

A. The zero-range force formalism 

For the purpose of obtaining the energy density functional in many-nucleon system, it is reasonable to begin with 
the nucleon-nucleon interaction of the following form: 

v * = E E, +EE E, , ,<•-,- ■ • • ■ ( 3 ) 

where the first and second terms of the right hand side are the two-body (two-nucleon) and three-body (three-nucleon) 
interactions, respectively. In this formalism all the interactions between each nucleon are summed up. With respect 
to the interaction used in density functional calculations, nucleon-nucleon interaction is usually considered up to the 
three-body term (cf. Skyrme interaction [To|). 

In this context special treatment (the short range approximation) is necessary to represent non-local potentials. 
First, the relative momentum between two nucleons, which is necessary to describe the interaction, is assumed to be 
represented by 

fe = |(v i -v i ), fc' = -i(v^-v;.), (4) 

where, again, i and j are the indices identifying the two nucleons. k' acts on the wave function to its left. For this 
representation, it is useful to remember the Einstcin-dc Broglic formula: p = hk under the quantum mechanical 
correspondence: p — > (fi,/i)V. This treatment is essential to obtain the zero-range interaction. 

The two-body interaction, which contains momentum dependence as well as spin-exchange contributions, is derived 
by assuming the polynomial expansion in powers of k and k 1 with the low-momentum cut-off up to quadratic terms 
(linear terms in k are forbidden by time-reversal invariance): 

Vi d {k, k') = i (l + x P a )S( ri - rj)) + |(1 + xxP^iSin - r 3 )k 2 + fc' 2 <5(n - + h(l + x 2 P a )k'5{r l - r 3 )k 
+ t -i\{S(a i ■ k')(aj ■ k') - (o-i • rT 3 )k' 2 }5{ ri - r,-) + 5{n - r i ){3( < 7 4 • k){rr 3 ■ k) - {rr t ■ a )k 2 }] 
+4f [3(<r. t • k')S(n - rj )(trj ■ k) - {a, ■ cr,)fc' ( 5(r. 1 - rj )k], 

(5) 

where to, ti, t%, t e and t a are parameters (following the traditional manner of representation), P a = ^(1 + ai ■ a 3 ) is 
the spin-exchange operator, and er.j are the Pauli matrices. The quadratic cut-off implies that the interaction shown 
in Eq (|5|) is valid to low-energy situations (cf. the fermi energy). The to term represents the central potential, the 
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t\ and t% terms stand for the non-local potential, and the t and t e terms are derived from the tensor interaction. It 
is important to introduce the spin-orbit interaction, which is an additional effect necessary to reproduce the nuclear 
structure. 

vff(k, k') = iW k' 6{n - r 3 ) (a. t + a 3 ) x k. (6) 
Another important ingredient is a two-body interaction with density dependence: 

vPP(k,k>) = ki+^p.m^)^-^). (?) 

This provides a treatment of medium effects including three-body and many-body forces. Note that the introduction 
of this force is due to the phenomenological reason. 

As a result, the nuclear interaction is represented by the zero-range Skyrme interaction: 

Vi,j(k, k') = v itj (k, k') + «£f (fc, fc') + v?f{k, k'), (8) 

where note that the included momenta k and k' are represented by Eq. (j4]). 



B. Hamiltonian density 



The Hamiltonian density is obtained by calculating the energy expectation value from the zero-range nuclear 
interaction jS]). The following points, which are illustrated in concrete discussion, should be noticed: 

• Fermionic relation is assumed in deriving the Hamiltonian density. 

• The boundary condition is necessary. 

• Differential operators appear in deriving non-local interaction part of the Hamiltonian density. 

In this section (Sec. IIIIB]) , wave functions are assumed to be sufficiently smooth. With respect to the many-nucleon 
system as a fermionic many-body system, fermionic relation (P r P a P T = — 1), which has already mentioned in Sec. lIII Al 
has to be taken into account: 

Vij(k,k') = ^Vi,j(k,k')(l - PrPaPr)- (9) 

Using Eqs. © and ©, the energy is obtained by calculating 

E = E to + E tl + E t2 + E t3 + E ti! + E to + Ew 

= ^Ei. m J ^(r^ m (r;-)K,(fc,fc')(l - P r PaPr)Mn)i>m(r i )dr i dr j dr , t dr' j , 



(10) 



where Et a , Et x • • ■ , Ew denote the energy arising from the terms with the coefficient to, t\ • • ■ , Wo, respectively 
In Eq. (|10p . ipi(ri)ip m (rj) plays a role of trial function including two single-particle wave functions. Furthermore 
we assume that there is no isospin mixing in the nuclear force. Since the expectation value E of the many-body 
Hamiltonian can also be written using Hamiltonian density H{r): 



E = J dr 3 H(r), (11) 



we obtain the representation for the Hamiltonian density. Among several terms, in order to sec some principal 
mathematical features (e.g., the appearance of differential operators) of the density functional formalism in nuclear 
physics, we examine the following energies, which are picked out from Eq. (1101) . 



E — E t „ + E tl + E t2 + E t3 , 



where this choice is also due to the pedagogical reasons. 
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Let ip be a single wave function, and introduce the densities. Before calculating the Hamiltonian density it is 
necessary to introduce the following densities. The nonlocal particle density and the nonlocal spin density are defined 
by 

S(r,r') = T,< T ,< T >, q Pq(r,r',<T,<r') < Cr '\ (T \ (T >= Ej.a.o-' ,q ^ ( r ' ' °' > «) CT ^ 0% °"> ?) 

using positions r and r', isospin q, and spins a and ct'. According to Engel et al. local densities are defined by 

p(r) = p(r,r')\ r=r >, 

T(r) = W-V'p(r,r')\ r=r ,, 

Sp(r) = s M (r, r')|r=r', 

J M (r) = -i(V M -V>(r,rO| P =r', 

^W = -i(V /1 -V;K(r,r')|r=w, 

T„(r) = V- V / « /1 ( r > r ')lr=*-'> 

where the indices /z and ^ indicate the Cartesian coordinate components x, y, and z. In particular pv in J M „(r) can 
take only unequal values to /z and v (cf. exterior product). The first two densities, that is, the particle density p{r) 
and the kinetic energy density t(t*) are scalar quantities, while the others, the spin density s /J (r), the current density 
jfi{r), the spin-current density J Mi ,(r), and the spin kinetic energy density T^fr) are vector quantities. The index q 
for isospin is used instead of r because it is used for the kinetic energy density. 

Galilean invariance restricts the form of the Hamiltonian density [2}. Galilean transformation replaces a single- 
particle wave function ip(r,a, q) by e lkr ijj(r, a, q). Thus the densities transform as 

p(r,r') e lk '( r ~ r '} p{r, r'), 
r(r) -> r(r) +2k ■ j(r) + k 2 p(r), 
j{r) -> j{r) + kp(r), 

because 

e~ ikr ^{r' ,a,q) e lkr ip(r,cr,q) = e lk{ - r ~ r '^{r' , a, q) tp(r,a,q), 

{Ver lkr 'i>{r', a, q)) (Ve lfcr ?/>(r, a, q)) 

= (-ike-' lkr ' '4>{r' ', a, q) + e- %kr ' V^(r', a, q)) (ike lkr ^(r, a, q) + e ikr Vij){r, a, q)) 
= feV fc ( r - p '>t/>(r', a, q)tp(r, a, q) + e ife (*--r') V^(r', a, q)Vip(r, a, q) 

-ike ik ^- r '^(r', a, q) (V^(r, a, q)) + ike lk ^ r '\V^{r' , a, q)) 4>{r, a, q), 

~(i/2){e- %kr '4>(r',a,q) (Ve ikr i)(r,a,q)) - (Ve" lfcr ' ip(r', a, q)) e lkr ip(r,a,q)} 
= -(i/2){e- lkr 'i>(r', a, q) (ike lkr ^(r, a, q) + e ikr V^jj{r, a, q)) 
-(-ike- ikr 'tp(r', a, q) + e - lkr ' Vtp(r', a, q)) e ikr ip(r, a, q)} 
= -(i/2){e lk ^- r ">4>(r', a, q)S7^(r, a, q) - e tk ^ r - r ">V^(r' , a, q)i>(r, a, q) + 2ike tk ^ r '^(r' , a, q) if)(r, a, q)}. 

The transformation leaves pr — j 2 invariant. Indeed, e ifc '( r ~ r 'p(r, r')\ r — r i = p(r) leads to 

p(r + 2k ■ j + k 2 p) - (j + kp) 2 = p(r + 2k ■ j + k 2 p) - j 2 - 2j ■ k ■ p - k 2 p 2 = pr - j 2 . 

Consequently j 2 cannot appear alone and the corresponding ingredient of the Hamiltonian density is pr — j 2 instead. 
We use the following substitution through integration by parts: 

(Vp(r)) 2 -p(r) (V 2 p(r)), (13) 
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which can be valid under a suitable boundary condition. In this section fSec. UHBj) . we assume that the Dirichlet-zero 
boundary condition is imposed. Note that nonzero additional boundary integral remains if we have no boundary 
condition. We have 

S7p(r,r')\ r=r > = Y,l, t7 ,qi ) l( r ^ a ^l)^^l-( r ^ a ^l)\r=r' = J2l.a,q M r , °, q)^J M r , ° , q) 

= £j,,r, g {i($ ; ( T 'i CT ' <l) V Mr, ^ l) + i>l(r, a, q)V^i(r, a, q)) 

+\({[} l {T,a,q)V'ipi{r,a,q) - ^i(r,a, q)Vrpi(r,a, q))} 
= ±Vp{r)+ij(r), 

V'p(r,r')\ r=r , = |Vp(r) -ij(r), (14) 

(V 2 + V' 2 )p(r, r')\ r=r , = J2t^ q Mr', <r, q)^Mr, <r, q) + Mr, a, g) V /2 ^(r' ', a, q)| r=w 
= Y.i.^ii.Mr', o, q)V 2 Mr, <r, q) + Mr, <r, q)V' 2 Mr' , <r, q) + 2V'$i(r' ) a, g)V^(r, a, ?)) 

-2VV/(r', a, q)Vipi (r, a, q)})\ r =r> 
= V 2 p(r) - 2r(r). 

For the historical milestones for obtaining the Hamiltonian density in many-nucleon systems, see Refs. [2L IT3|. 



i. Force depending on the parameter to 
First, for the term with the coefficient to, we have 

(l + X Pa)(l-PrPaPr) 
= (l+X P<r)(l-PaP T ) 
= 1 + X P a - (x P 2 + P a )P T 

= 1 + ±X - |(1 + 2x )P T + \xQ<Ti<Tj - \<Ji<JjP T , 

where we can take P r = 1 due to the fact that the zero-range force acts only in S-waves. It leads to the energy 

Et = ~2 Yjl,m 

J Mr'i)M{r'j) {l + \x - |(1 + 2x )P T + ^a.cTj - \cTiajPr} MrM'mir J )dr l dr ■ j dr' i dr' j \ ri=rj=r -- r '. 
= Jilt (1 + k x o)p(n,ri)p(ri,r() - f (| + x )5 quq2 p(r i ,r^p(r i ,r^ (15) 
+i%XQs(r i ,r' i )s(r l ,r' l )-if5 qim s{r l ,r[)s(r l ,r' i )}dr l dr dr'^r'- \ Ti=Tj = r ,. 

= / dr 3 {f (1 + \x Q ) P {rf - f (i + aro) E 9 P^) 2 + ^o*(r) 2 - | E 9 ^(r) 2 } , 

where P r is reduced to 5 qiy q 2 by assuming that there is no isospin mixing in nuclear force (qi identifies the isospin of the 
single particle state i). The specific representation shown in the second line of Eq. ([15)) , which leads to the derivation 
of the Hamiltonian density, has been used in physics, where the property of ^-function: J S(a — r)f(r)dr = /(a) and 
therefore 



S(r - r')f(r')dr'dr = J f(r)dr 



is taken into account. 



2. Force depending on the parameter ti 



Second, for the term with the coefficient ti, we have 
(f + x 1 P a )(k 2 + k' 2 )(l-P r P a P T ) 

= {(U^i v,)) 2 + (-^(v< - v;-)) 2 }(f + x 1 p„){\ p a p T ) 

= {-|(v? - v i v j - v,v, + vf) - i(vf - vjvj - v;.v^ + vf )}(i - Xi p 2 p t + Xl p a - p a p T ) 

= -i(V 2 + V 2 + Vf + Vf - 2ViV 3 - - 2V' 2 V;-){f - Xl P T + \{ Xl - P T ){\ + o^)} 
= -i(V 2 + V 2 + Vf + Vf - 2ViV,- - 2V' 2 V;-){f + \xi - {x x + |)P T + - Pr)^-}, 
where we can take P r = 1 due to the fact that the zero-range force acts only in S-waves. It leads to the energ 

E tl = -&Zt, m J MrW m (r>)(V 2 + V 2 + Vf + Vf - 2V.V, - 2V{V$) 

{1 + \x\ - (xi + \)P T + \(xi - P T )a i a j }il>i(r l )i) m {r J )dr i dr J dr[dr' ] \ r . %=r ..= rli = rli 

= -H /(V 2 + V? + Vf + Vf - 2V i V ? - 2VJV$) 

{(I + \x x )pir i ,r^p{rjy j ) - (x x + \)5 quq2 p{r u r^p^j, r'j) + \xxs{ri, r'^s{rj, r\) 
-^qi,q2 s ( r i, r i) s ( r 3, r j )}dr t drj dr[ dr^ \ rz=rj =r , =r , 

= -£Jdr 3 { (1 + ±xi){[(V 2 + Vf)p(n,rOMr j7 r;) + p(n, r ? ')(V 2 + Vf)p( rj ,r>) 
-2V 4 p(r l ,rOV J p(r„r;.) - 2V' i p(r i ,r! i )V' j p(r j ,r' j )} 

~(xi + ^)EJ(V 2 + VfKKrOWr,,^.) + p 9 (r 4 ,rQ(V 2 + V*)p,(r 3 -,r<) 
-2V i p,(n,rOV jP? (r J ,r;.) - 2V<p g (n, r^V^r,, r<)] 
+Ix 1 [(Vf +Vf) S (r ! ;,r,:)} S (r J ,^.) + S (n,r ? ')(V 2 + Vf) S (r„r;.) 
-2V iS (n,r<)V iS (rv,r;.) - 2Vj*(r if rOV>(r 3 -, r<)] 

= -ft / ^ 3 { (1 + ^i)[2(V 2 p(r) - 2r(r))p(r) - 2(§Vp(r) + *.?(r)) 2 - 2(|Vp(r) - ij(r)) 2 ] 
-(xi + \) Ej2(V 2 p 9 (r) - 2r q (r))p q (r) - 2(±Vp ? (r) + i j (? (r)) 2 - 2(iVp 9 (r) - ij q {r)) 2 ] 
+ixi[2(V 2 s(r) - 2T(r))s(r) - 2(±Vs(r) + U{r)] 2 - 2(±Vs(r) - *J(r)) 2 ] 
"3 E 9 [2(V 2 S ,(r) - 2T q (r))s q (r) - 2{\Vs q {r) + i J q (r)) 2 - 2(±Vs 9 (r) - iJ q {r)) 2 ]} 

= Idr 3 {^[-(1 + ix 1 )p(r)V 2 p(r) + ( Xl + §) £ g P 9 (r)V 2 p 9 (r)] 

+£[(1 + \ Xl ){p{r)r{r)-j{r) 2 ) - (xi + \) £ g (p 9 (r)r,(r) - j q (r) 2 )} 

-^[ XlS (r)\7 2 s(r)^j: q s q (r)V 2 s q (r)} 

+ t i[x 1 (s(r) • T(r) - J 2 (r)) - £,(s,(r) • T,(r) - J 2 (r))]}, 

where we use the relations shown in Eq. ()14|) and the reduction: P T — > S qi m . 



3. Force depending on the parameter ti 

Third, for the term with the coefficient t 2 , we have 
(I + x 2 P a ){k ■ k')(l - P r P a P T ) 
= i-Ti^i - Vi) x ^(V^ - W' j )}(l + x 2 P rr )(l+P lT P T ) 

= i(ViV< - Viv;. - v,-v< + v, v;.)(i + x 2 p T + X2 p a + p a p T ) 

= i(v,v^ - Viv;. - v 3 v< + v,v;.)(i + s 2 p T + |(x 2 + p T )(i + v^)) 

- KV<V< - V,V;. - VjVJ + V,V;.)(1 + ^-2 + (x 2 + \)P T + \{x 2 + Prfatrj), 
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where we can take P r = — 1 since odd powers of fe are treated. It leads to the energy 



(1 + \x 2 + (x 2 + \)P T + 5(^2 + Pr^iVj^iird^mir^dridrjlr.^.^ 

= f /(v*v{ - v,v;. - v, v 2 + v, v;.) 

{(1 + \x 2 )p{r i ,r$p{r j y j ) + (x 2 + 5)V,g 2 P( I 'i,^)K r j 1 r j) 
+iaj 2 s(r i ,r t Os(r J -,r^) + \b~ quq2 s{r u r<)s(r,-, r<)} dr t dr 3 | ri=rj=r j =I .- 

= | fdr 3 [(l + ±z 2 ){p(r)r(r) - (|Vp(r) - i(r))(±V P (r) + j(r))} 
+ (^2 + I) E q {P q (r)r q (r) - (§Vp,(r) - j,(r))(|Vp,(r) + i,(r))} 
+ i.T 2 { S (r) ■ T(r) - (iVs(r) - J(r))(|V • s(r) + J(r))} 
+ 5 E, KM • r,(r) - (^Va,(r) - J,(r))(|V • «,(r) + J q (r))}] 

= %J dr 3 [(l + ix 2 ){4p(r)r(r) + p(r)V 2 p(r) - 4j(r) 2 } 
+(*2 + |) £,{4p,(r)r g (r) + p 9 (r)V 2 p(r) - 4j g (r) 2 } 
+.T 2 {2s(r) • r(r) + §a(r) • V 2 s(r) - 2J(r) 2 } 
+ E 9 {2s ? (r) • T,(r) + \s q {r) ■ V%(r) - 2J(r) 2 }], 



where we use the relations shown in Eq. (fT4"]) and the reduction: P r — > 5, 



Fourth, for the term with the coefficient t 3 (the density-dependent term) , we have 



where we can take P r = 1 due to the fact that the zero-range force acts only in S-waves. It leads to the energy 



E ts = ||E[,m / M r 'i)$m(rj) {1 + ix 3 - i(l + 2X 3 )P T + ^X^fTj - \(T l (T : jP T } 

+^x 3 s(r l ,r , l )s(r l ,r! i ) - ^5 quq2 s{r il r[)s{r l ,r[)p{^^) a }dr l dr j \r l =r J =r[=r' ] 
= /{&(! + ^3)p(r) 2+Q - H(i + a*) E 9 P q (r) 2 p q (r) a + ^s{rf p(r) Q - f| E g *,(r) 2 p,(r)"} dr 3 , 



^. Force depending on the parameter £3 



(1 + z 3 P CT )p(E£p)«(l - P r P CT P T ) 

= (1 + x- 3 p (T )(i - p ff p T )p(n+Si)«» 

= {1 + x 3 P a - P a P T - x 3 P*P T }p( 



{1 + \x 3 - |(1 + 2x 3 )P T + lasffiffj - \(r l a 3 P T }p{^p-) 



where we use the reduction: P T — > S, 



5. Hamiltonian density of interacting many-nucleon systems 
The expectation value E of the many-body Hamiltonian is written by an integral of the Hamiltonian density H(r). 

E = E t0 + E tl + E t2 +E t3 = [ dr 3 H(r). 
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The Hamiltonian density is obtained by collecting the results: 

H <r) = \[&J + W + \^)p 2 ~ f (I + *o) E g p\ + 1 ^s 2 - 1 E 9 * 2 

-%(1 + ±Xi)pAp + ^.(xi + i)E,P 9 Ap 9 + ^(1 + \x{]{pr-f) - \{x x + \)Y. q iP q r q - P q ) 
-^ XlS As + %Y, q s q ^s q + ^ Xl {s-T-^)-^Y, q {s q -T q ^J 2 q )] 

(16) 

+|f(l + ±z 2 )pAp + |(1 + ix 2 )(pr - j 2 ) + £,{f|(a: a + i)p 9 Ap + t -f(x 2 + i)(p ? r g - j 2 )} 
+||x 2 s • As + f z 2 (s • T - J 2 ) + EJi^g ■ As 9 + |(s, • T q - J 2 )} 

+ \^)p 2+a - ff (I + »a) E, ^p? + li^V - & E, « 2 p1 ; 

it is half of the standard Hamiltonian density 0, fl3j , which can be understood by applying the variational principle 
to the first term of the right hand side (the contribution from the kinetic energy), where 1/2 for the interaction 
part appears to avoid calculating the interaction twice. The justification of 1/2 (in the term arising from the kinetic 
energy) is illustrated in the next section (Sec. IIII C[) . It is notable that differential operators appear due to the 
non-local interaction. Interacting many-nucleon systems cannot be described well if differential operators are not 
utilized (cf. local density approximation). 



C. Energy density functional 

The components in the effective interaction are obtained based on the variational principle. This corresponds to 
the procedure of obtaining the functional representation of an effective Hamiltonian. The following points, which are 
illustrated in concrete discussion, should be noticed: 

• Once the energy is given, it is not necessary to have the fermionic relation or the Slater-determinant formalism 
to derive the effective interaction (the effective Hamiltonian). 

• The boundary condition is indispensable. 

• Regularity of wave functions is required; i.e., L 2 -space is not necessarily sufficient. 

• For mathematically strict treatment, the variational principle should be considered not only in the real 
Hilbcrt/Banach spaces (by varying "real" densities) but in the full complex Hilbcrt/Banach spaces (by varying 
"complex" wave functions). 

Since the application of the variational principle includes differentiation, more careful treatment is required (differen- 
tiability and so on). Furthermore, the application of variational principle means a process of obtaining a functional 
representation in some complex functional spaces only from the real function (i.e., energy). In order to present the 
essential treatment of the variational principle, the discussion is developed for a spin-saturated and charge-conjugate 
nucleus (i.e., p/2 = p p = p n ). In this situation Eq. (|16[) reduces to a model Hamiltonian density: 

H(r) = fir + t 2P 2 + hp* + U(pr - j 2 ) + t 5P Ap, (17) 

where a is a positive rational number, and {ti} is a parameter set related to the parameter sets {ti} and {x{} 0]. 
Here, based on this model Hamiltonian, generalized higher order terms and terms with fractional powers are also 
treated, where we need to discuss within suitable functional spaces for a strict treatment of the variational principle. 
First, the term arising from the kinetic energy (the homogeneous term of Schrodinger equations) is treated in Ex. 
IIII C II Next, the derivation of terms arising from the interaction energy is demonstrated. 

1. Linear and nonlinear Laplacians 

Let Z/ P (f2) be the space of functions on f2 which are LP for the Lebesgue measure, where L p {Vi) are not (infinite- 
dimensional) Hilbcrt spaces but (infinite-dimensional) Banach spaces if p is not equal to 2. Let p be an even number 
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satisfying p > 2 and ip be included in VT~ 1,p (f2) with a open bounded set il C R 3 , where W 1,p {Vt) is the space of 
functions in L p {Vl) whose distribution derivatives of order < 1 are in L p (f2) (for the "mathematical" distribution 
refer to the textbook of functional analysis) . For our purposes it is necessary to introduce such a generalized differential 
in association with the spaces of integrable functions. Note here that the 5-function, which cannot be defined using 
usual functions (C^-function and so on), is successfully defined in the sense of a distribution. In this sense the 
zero-range interaction formalism is based on the theory of distribution. 

First of all it is useful to consider a simple case when the energy is represented by 

£(xjj) = - I dr 3 \V^\ p , 



PJn 

where the corresponding Hamiltonian energy density is H(ip) = -|V-0| P . The differential operator is considered in the 
sense of distribution; this is true to other cases shown in this section (Sec IIII Cj) . By minimizing this energy (i.e., the 
Gateaux differential: £'(V>) = A is considered), we obtain 

A{i>) = -v-(|W| p ~ 2 W0- 

Indeed, for tp and </> in W 1,p (Vl) : 

|W> + \<j)\P - \Vyj\P = /„* £|V^ + tAV0pd* 
= /of(l v ^ + *AV0| 2 )P/ 2 di 

= Jo f dW' + UV0| 2 )(P- 2 )/ 2 f |VV> + t\V<j)\ 2 dt (18) 
= Jo f (I W' + aV0| 2 )( p " 2 )/ 2 ■ 2Re{(V?/> + tAV0) ■ AV0}dt 
= p\ J,, 1 \Vip + t\V(j)\P- 2 ■ Rc{(V^ + tAV0) ■ V4>}dt 

is true. It follows that 

£(v+A0)-gw = dr 3 j q i | V ^ + av0|p- 2 Rc{(v^ + t\V(f>) ■ V0}dt. 

Due to A (A > 0), 

|V?/> + aV0| p - 2 Rc{(W> + tAV0) ■ V0} -> |V?M P " 2 Re(VV' • V0) a.e. 

is valid ("a.e." means "almost everywhere", whose mathematical definition should be referred to the textbook of 
functional analysis. Furthermore, due to A — > 0, the following inequality is valid: 

\\Vip + t\V(t>\ p ~ 2 Re{(V^ + t\\7<p) ■ V0}| < |V^ + tAV0| p_1 |V0| < (|V^| + |A||V^|) p - 1 |V^|, 

where it is worth noting here that the last term is integrable. We have 

±S(tp + A0)| A=Q = lim A ^ £ ^+^)- £ W = J n dr s |V7/>|f- 2 Rc(V^ • V^) = Ref Q dr 3 |V^| p - 2 V^ • V^. 

If ip and 4> are sufficiently smooth, 

f |VV>r _2 VV • dr 3 = f |W| p - 2 ^-0 dS+ [ A(ip)4> dr 3 
Jn Jan "l* Jn 

follows from the integration by parts, where dfdy, means the outward normal differential operator. Although we 
have discussed the limit-process in a roundabout sort of way, it is based on the fact that liniA->o £ ^ +x< ^~ £ ^ = 
— j n dr 3 \\7ip\ p ~ 2 VipV(f> cannot follow from some simple treatments, because the left hand side is real-valued and the 
the right hand side is complex-valued. It means that the differentiability of the energy cannot follow naively (for the 
differentiability, see Appendix (Sec. EH). This fact is true for all the cases treated in Sec. IIII Cl 

The final treatment actually depends on the boundary condition. When the Dirichlet boundary condition is adopted 
(for the homogeneous part, Schrodinger operator is considered), Vy 1,p (f2) is replaced by W ' p (O), and we have 

IWP^W • V0 dr 3 = [ A{ip)4> dr 3 = (A(ip), <f>). (19) 
Jn 
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That is, for / G W~^ p ' {Vt) = W 1,p (n)* and ip £ Wo' p (fl), the equality: f Q |V^| P_2 V^ ■ V0 dr 3 = (/, 0) is valid to 
any ip G Wg' p (f2), which is an generalized solution of 

f .A(vo = / in n, 
I V = o on an. 

Meanwhile when the Neumann boundary condition is adopted, we have Eq. (fTi?|) . However the detail is different from 
the case with the Dirichlct boundary condition; for / G W l ' p (Cl)*, the equality: J n \\7ip\ p ~ 2 Vip ■ V<^> dr 3 = (/,</>) is 
valid to any ip G W 1,P (Q), which is an generalized solution of 

A(tl>) = f inn, 

U = o onan. 

When the periodic boundary condition is adopted for cuboid n, W 1,P (Q) is replaced by W^(Vl) 1 which is the space 
of restrictions to n of periodic functions. We have Eq. (fT9|) using the similar argument. As a result, 

±£(ip + A0)| A=O = Re j n dr 3 \Wtp\ p - 2 Vip ■ V0 = Re J n A(ip)4> dr 3 = Re(A(tp) , <j>) (20) 

is valid to the Dirichlct, Neumann and periodic boundary conditions. The operator A with p > 2 is called p— Laplacian 
or nonlinear Laplacian. If we take p — 2, the relation between A, H and £ is that between the Laplacian, the 
corresponding Hamiltonian density and the corresponding energy. The case with p = 2 has the direct connection to 
the standard kinetic energy term. Note that, for this discussion, it is sufficient that ip is included in W 1,p (n), and it 
is not necessary for "0 to be a component of the Slater determinant. 

In case of many-particle situations there is no need to introduce extraordinary treatment. Let the many-particle 
wave function be 

TV 

i=i 

It is trivial that the energy becomes 

N 

£nW = - / dr 3 VIV^p, 
pJn 

where V* acts on ipi, and the term in the effective Hamiltonian 

N 

i=i 

follows. Here the procedure of obtaining a term in the effective Hamiltonian from a given energy consists only of 
Gateaux differential, so that the similar procedure is also valid if we assume many-particle situations. When we 
consider the linear situation p = 2, the term is equal to X)^Li( — ^iV'Oi corresponding to the homogeneous term of 
the typical type of many-body Schrodinger equations. In the following we do not discuss the many-particle situations 
explicitly, but ,as is the previously shown, they trivially follow once the corresponding one-particle situations are well 
understood. 



2. Nonlinear interactions depending on the parameters ti and ti 



Again let p G Z (Z: a set of all integers), be an even number satisfying p > 2. We consider the energy represented 
by 

£dA1>) = - f dr3 W p > ( 21 ) 



V Jn 

where the corresponding Hamiltonian energy density is Hd,z{^) = By minimizing this energy (i.e., the 

Gateaux differential is considered), we obtain 
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where this is one of the most typical ingredient of nonlinear Schrodinger equations. Indeed, for ip and <fi in L P (Q), 

\ip + A0|p - = Jo §- t |V + t\<j>\Pdt = P X Jo IV- + U0|f- 2 • Re{(^ + t\<j>) ■ 4>}dt (22) 

is true. 

\i> + t\(f>\ p - 2 Rc{(i> + t\<p) ■ $) -> |V>r -2 Re(V> • 4) a.e. (23) 
is valid due to A— > (A > 0). In the same manner, we have 

&£ D ,zW + W)\ x=0 = linu^o e °.*W+W-SoM1>) = /o dr 3 |^ r 2 Rc(V ,.^ ) = Rc ; odr 3 
As a result, we have 

££ D ,z(4> + A0)| A=O = Re J n dr 3 • 4> = Re/ n ^(V)^ dr 3 = Re(.F AZ (V>), # (24) 

Nonlinear Schrodinger equations only including the interaction Tu,zi^) with p = 4 are known as one of the most 
typical type of nonlinear Schrodinger equations (cf. Ginzburg-Landau formalism). 

Let us consider the energy £d,q when p £ Z in Eq. ([21]) is replaced by q € Q (Q: a set of all rational numbers). 
Nonlinear interaction with q £ Q, which is simply called density-dependent force, is suggested to be indispensable to 
explain experimental results: 

where q satisfies q > 2. For example, q satisfying "(g — 2)/2 = 1 + 1/6" has been proposed as a possible candidate [l|, 
as well as that satisfying (q — 2)/2 = 1 + 1/4 0. Note here that the "density-dependent" force is a technical term 
for the force J 7 d j q(V') satisfying q ^ 4. If we simply use the previous discussion shown in Eqs. (|2"2")l and (f2"3")l . it is 
necessary to introduce the L 9 -space of the fractional power. Here is not problem to define L 9 -space with noninteger 
and such L 9 -space satisfying 1 < q < oo holds the property of Banach spaces. Accordingly the discussion shown in 
Eqs. (f2"2"|) and (j2"3")l are valid even in this case, and we have 

&£d,q(1> + A0)| a=o = Re/ Q dr 3 W~ 2 ^ • - Re/ n T d ,qW)4> dr 3 = Re(F D>Q (i>),4>). (25) 



3. Nonlinear interaction depending on the parameter i± 

Again, let p G Z (Z: a set of all integers), be an even number satisfying p > 2. We consider the energy represented 
by 

mw = - / <fr 3 iwr 2 M 2 , 



where the corresponding Hamiltonian energy density is Hk,z(iP) = -|V'0| p ~ 2 |V-'| 2 - By minimizing this energy (i.e., 
the Gateaux differential is considered), we obtain 

tk.z{^) = - • awMvfwo + |iv^r 2 ^. 

Indeed, for ip and <j> in W 1,p (f2), 

|V(V> + A0)|f" 2 |V + A0| 2 - |W|p- 2 |i/;| 2 = J,, 1 f {|V(^ + tA0)|^ 2 |^ + U0| 2 } A 
= A [(p-2)|V(?A + a0)|P- 4 |V + iA(/)| 2 • Re{V(V> + £A</>) • V^} (26) 
+2|V(^ + U0)|p- 2 • Re{(^ + *A</») • 0}] 

is true. 

(p- 2)|V(^' + £A0)| p - 4 |V + *A0| 2 •Rc{V(^ + U0) • V^} -> (p - 2)| Vijj\ p - 4 \ip\ 2 Re{Vi< ■ V0) a.e., 
2\\7(i' + t\(t))\P- 2 • Re{(V> + *A<£) • 0} -> 2|VV>| p_2 Re(V' • <£) a.e. 

are valid due to A — >• (A > 0). We have 

jfc&r.zty + A0)| A=Q = Re± / n dr 3 {(p - 2) | W |^" 4 1^ | 2 ( • V0) + 2|VV>r 2 (V> • 0)}- 



(27) 
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If ip and <fi are sufficiently smooth, 

/ |VV>r" 4 |V| 2 VV" Vj> dr 3 = [ |VV>r~ 4 |Vf 2^<MS--Re / V-(|VV-r" 4 |V| 2 V^)-0dr 3 . 
Jn Jon "fx J n 

As a result, by taking into account the boundary condition: the boundary conditions shown in Sec. 1111 C 1| we have 
A^ zW , + A0)| A=o = -Re/ n dr 3 ■ (| WMVfVV) • + Re/ dr 3 ||V^r^-0 

= Re f n FK,zU>)$ dr3 = M^K,z{i>)A)- 
Note that this interaction including differential operators cannot be bounded on L 2 (S7). If p = 4, J^.z is reduced to 

^,zW= -iv-dV'I'v^ + ilv^i 2 ^. 

Similar to Sec. 1111 C 21 nonlinear interaction with q £ Q is derived from the energy Ek,q, which is obtained by 
replacing p G Z in -Er-,z by q G Q. 

^,qW) = - • (|vv>|*- 4 M 2 vvo + \ |W>r 2 v>. 

where g satisfies g > 2. If we simply use the previous discussion, it is necessary to introduce the W 1,9 -space of the 
fractional power. Using the interpolation of Sobolev spaces, there is no problem to define W 1,9 -space with noninteger 
q, and such W^'-space satisfying 1 < q < oo holds the property of Banach spaces. Accordingly the discussion shown 
in Eqs. |26|) and (|27|) are valid even in this case, and we have 

&Sk,q(iI> + A0)| A=0 = Re / n dr 3 \V^~ 2 ^ ■ = Re /„ T K:Q (^ dr 3 = Rc(T K:Q (^), 0). (29) 



4- Nonlinear interaction depending on the parameter t^ 

The nonlinear term represented using the Laplacian operator is considered. We consider the energy represented by 

£ L {i>) = \Jjr 3 \i>\ 2 ^\\ 

where the corresponding Hamiltonian energy density is rlhi^) = \ IV'PAIV'I 2 - By minimizing this energy (i.e., the 
Gateaux differential is considered), we obtain 

T L {$) = (AH 2 >/>. 

Indeed, for ip and 4> m W 2 4 (Q), 

\ij> + A0| 2 A|^ + A</>| 2 - H 2 A|Vf = /„ §- t {\ip + t\<j>\ 2 A\ip + U0| 2 } dt 
= 2A [(A|V> + U0| 2 ) Rc((^ + tA0) ■ <j>)+ReA((ip + t\tj>) ■ <p)\ip + t\<j)\ 2 ]dt 

is true, where 

(A^ + tA^ 2 ) Re((V> + iA</>) ■ 4>) -> (A|^| 2 ) Re(-0 • <j>) a.e., 

A(|V + tA^| 0)|V> + ^| 2 -> A(|V#) IV-I 2 a.e. 

are valid due to A — ► (A > 0). We have 

^£ L (V + A0)| A=Q - | / n dr 3 {(AH 2 ) Rc(^ ■ 0) + RcA(-0 ■ 4>) |^| 2 }. 

If ip and ^ are sufficiently smooth, by taking into account the boundary condition: the boundary conditions shown in 
Sec. 1111 (J II we have 

/ n A(M </>)|?/>| 2 dr 3 = -/ n V(|V| 0)V|V| 2 dr 3 = J n A|V>| 2 Re(V 0) ^ 3 . 

Consequently, 

^£ L (V + A0)| A=Q = Re j n dr 3 {(A|V>| 2 ) ^ • 0} = Re ! n F L {^)4> dr 3 = Re(Ji(V)^). (30) 
Note that this interaction including differential operators cannot be bounded on L 2 (il). 
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5. Nonlinear interaction arising from current (depending on the parameter ti) 

The momentum density j , which is defined by j = — | (ipVip — tpWip), plays a significant role more than one physical 
quantity to be utilized to describe the interaction. For a master equation: id t tp = (1/2) (—A + V)ip in a Hilbcrt space 
L 2 (ft), 

= §(i(-A + V)tp)ip + |^(-i(-A + V)ip 

= f ((-AV# - <K-AV)) + - VW))- 

If V is a self-adjoint operator in L 2 (Vt), then J^(V>V) = (^((-AVOV - V>(-AV0) = («/ 2 ) v • ($W ~ ^V$) follows, 
and 

-(^) + V-j = (31) 

is obtained. This equation, which means the conservation of total particle density (represented by J dr 3 \ip\ 2 ), is known 
as the continuity equation. That is, the momentum density j plays a role of current. 
We consider the energy represented by 

£j(i>)= [ dr 3 j 2 = ^f- f dr 3 (^V^-^) 2 , 

where the corresponding Hamiltonian energy density is T-Lj{tp) = —\($>Wij) — ^V^) 2 . By minimizing this energy (i.e., 
the Gateaux differential is considered), we obtain 

FjW) = -2*{2(V-j> + i-W. 

Indeed, for -0 and (f> in W 1,4 (0), 

{(tp + X^)V(ip + \4>) - (?A + \<p)V{ip + X(j))} 2 - ($vv> - ^V$) 2 
= f {(V' + tA0)V(V' + iA</>) - (-0 + t\<l>)V(il> + t\<l>)} 2 dt 
= Jo d£ [2{(fTi^)V(^ + f A0) - + U0)V(^ + U0)} 

&{ty> + tA0)V(v> + tX4>) - (V- + tx<p)v(^ + tx<f>)}] 

= 2X Jq 1 di [{(VH r W)V(V' + a^) - (tp + tX<j))V{i> +1X4)} 
{<j)V(^ + tX<f>) + (V' + tX4>)W<f) - <t>V(ip + tXcf>) - (ip + tXcj))V4>}] . 

It follows that 

iA±kM=£lM = _ i J n dr 3 ji ^ [{(^ + tA0)V(^ + iA0) - (V- + U0)V(^TlA0)} 

{^v(?a + tx<t>) + (ip + t\<t>)v<t> - <jN{tp + tx<j)) - a< + tx<j))y^}\ . 

Due to A -> (A > 0), 



-{{il} + tX(j))V{i) + tX(j)) - {ip + tX(/))V(ip + tX(f>)} ->■ - {^VV» - ^V^} = 2ij a.e. 



{<£V(V> + iA0) + OA + tA0)V0 - 0V(t/> + <A0) - + iA0)V<?i} -4 0V?/> + - 0V?/> - V'V^ a.e. 
are valid. If ip and <J) are sufficiently smooth, we have 

&£j(1> + A0)| A=Q = - /„ dr 3 (ij){(V^ + ^ - (V</0<£ - 
= ~ Jan ' dS + fan { U ' v )$<t> dS 

7iJA6 + ii(X7iJA6 - i(X7 ■ i)ib6 - ii(X7i/Ad>\ 

(32) 



+ / n dr 3 {-ij{V^)4> + i(V • i)V^> + ij(yi>)4> + U'(VV)^ - »'(V • j)M - «j(W>)0} 



~ Jan ' d5 + J an i(j • v)if)(f> dS 
+ f n dr 3 {-2ij(Vip)(j> + 2ij{Vip)(j) + i(V • j)$<j) - i(V • j)i/)<j>} 

- f 9Q i(j ■ v)n dS + j m i(j ■ v)W dS + j n dr 3 {-4Rc^(VV)0) - 2Rc(i(V • j)^)}, 
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where v means the outward normal vector for the dfl. By taking into account the boundary condition: the boundary 
conditions shown in Sec. IIII C 1[ the integral on the boundary surface cancels. Therefore 

±£j{ij> + A0)| A=Q = Re J Q dr 3 [-2*{2j(W)0 + (V • = ( Jj (</>), 0), 

where it is very important to note that the second term of the right hand side of the following equation: 

(^W,^) =Re / -{j(Vvj)0+(V-j)^}dr 3 -Re [ -(V-j)^dr 3 (33) 

is missing in the standard formalism shown in the Appendix B of Ref. Q . This term arises from the careful treatment 
of the integration by parts (see Eg. ([32]) ). Note that this interaction including differential operators cannot be bounded 
on L?(VL). 



6. Effective Hamiltonian 

The corresponding effective Hamiltonian for the Hamiltonian density shown in Eq. (|17[) is obtained as 
H e s(r) = -ii2A + t 2 Ap + ? 3 2(1 + a)p a 

+i 4 {-2V • (h/fVVO + 2|VV|V + 2i{2j(W) + (V • + t 5 4Ap 

(34) 

= 2[-ti A + 2t 2 p +{2 + &)t 3 p 1+a 

+< 4 {-V • (|V| 2 V^) + |WfV + »{2j(W) + (V • .?>}} + 2t 5 Ap], 

where we should pay attention to the coefficients with signs. Two points should be noticed: as is seen in t\, if we 
obtain the coefficient in the effective Hamiltonian as — h 2 /2m, the corresponding coefficient t\ in the Hamiltonian 
density should be h 2 /4m (it requires the modification of the standard choice of coefficients contained in the Hamil- 
tonian density 0, HH); as is already pointed out, it is necessary to subtract 2i(V • j)tp (compared to the standard 
parametrization shown in the Appendix B of Ref. (H) in order to have a complete set of terms arising from j 2 in the 
Hamiltonian density. For the former point, it is useful to remember that 

/ rdr 3 = / -(A^dr 3 ; T =\Vip\ 2 
Jo. Jn 

does not follow from the variational principle, but from the integration by parts with a suitable boundary condition. 
Although the representation of the standard Hamiltonian density is not correct 0, [l3| , the standard effective Hamil- 
tonian based on such a Hamiltonian density is exactly the same as the effective Hamiltonian obtained here (except 
for the terms arising from j 2 ). For the latter point, the Galilean invariancc has to be broken if the missing term is 
not included (for the Galilean invariance, see also the corresponding discussion in Sec. IIII B|) . 

In this section the discussion has been developed in LP spaces, where p is not necessarily equal to 2. Indeed, 
p (y^ 2) is necessary to consider nonlinear problems; e.g., for tp £ L 2 (fl), \ip\ 2 is included in L 1 (f2) at least, and 
therefore we are not sure if \ip\ 2 ip is included even in L 1 (J7). Note again that L p (fl) are not Hubert spaces, if p is not 
equal to 2. In such situations (p ^ 2), the inner product is not equipped, so that we cannot discuss the orthogonality 
of functions and so on. Here is a difficulty of considering nonlinear problems. 



IV. REMARKS ON THE ADDITIONAL FORCES 



A. The Coulomb force 



Only protons hold charge, so that they interact by the Coulomb force. Charge is assumed to be equal to the 
probability distribution of proton (cf. form factor). The Coulomb force consists of the direct and the exchange parts. 
The direct part of the Coulomb energy is represented by 

Et -iff 
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Fictitious empty space 




FIG. 2: An actual numerical calculation of the Coulomb interaction is presented in the coordinate space when the periodic 
boundary condition is imposed. A white-coloured square means the computational cell. The solution is obtained by two fast 
Fourier transform operations in the enlarged region with periodic boundary condition. This treatment realizes the isolated 
charge distribution in the computational cell. 



using the proton density, where the corresponding Hamiltonian density is equal to ^- J dr 3 | ^'_^ | ■ On the other 
hand, the exchange part of the Coulomb energy is approximated by means of the Slater approximation [TlT ] as 

where the corresponding Hamiltonian density is equal to — (3e 2 /4)(3/-7r) 1/ ' 3 /9p(r i ) 4 / 3 . Note that the exchange part of 
the Coulomb force has the similar form as the term with the coefficient £3, so that the treatment of obtaining the 
corresponding part of the effective Hamiltonian is similar to the cases with E^ 3 . 

Let us move on to the numerical calculation of the Coulomb force. Because of the long-range property of 
the Coulomb interaction, the periodic boundary condition is not necessarily appropriate, but the potential has 
to go to zero at infinity ("isolated charge distribution"). In practice this case is solved usually in one of two 
ways, either by obtaining boundary values from a multi-pole expansion or by Fourier techniques embedding the 
computational cell in one large of empty cells (Fig. 2). In the periodic case being useful for astrophysical situation, the 
jellium approximation is used corresponding to a constant background density of electrons cancelling the total charge. 



B. Pairing force 

The pairing interaction is an important ingredient of quantum many-body systems. This interaction combines two 
fermions into one boson, so that condensation can take place as a new feature. For example the superconductivity 
follows from the pairing interaction. 

In order to introduce the pairing interaction to the density functional theory, it is necessary to have a functional 
representation of the pairing-interaction field. There have been proposed two methods of introducing the pairing 
interaction in many-nucleon systems @; one is the HFB (Hartrcc-Fock-Bogoliubov) approach, and the other is the 
BCS (Bardeen, Cooper, and Schrieffer) approach. 

Following the review article [l2j], here we introduce two pairing- field of BCS type: 



<q 4 

and 

pair 



E pai r=J2 q V i-fdr 3 Xq (r) 2 (35) 
E q V -T fdr 3 [{l-(^?) 7 }x,W 2 ]. (36) 
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where 

Xq(r) = ^ WgU & v & \iJj & (r)\ 2 

denotes the pairing density with the phase-space weight w&, occupation amplitude v&, and non-occupation 
amplitude u& = \/l — V | . In Eqs. (|55|) and (j5o]) . u g and uo. g are strength parameters, and po is the nuclear 
saturation density, typically p — 0.16 fm 3 . In contrast to Eq. (pjj) . Eq. d3U) additionally includes the density 
dependence. The surface profile of the pairing interaction is controlled by the corresponding parameter 7, whose 
standard value is equal to 1. Note that the strong pairing takes place near the nuclear surface, so that Eq. ([55]) is 
expected to describe the pairing field better than Eq. (|35p . The corresponding Hamiltonian densities for Eq. ([55)1 
and Eq. (j56"|) are equal to J2 q ^tXq{ r Y ancl J2 q li 1 - ( £ ]T~) 1 }Xq( r ) 2 }, respectively. As is readily seen, the 
treatment of obtaining the corresponding parts of the effective Hamiltonian are similar to the cases with Ef 2 and Ef 3 . 



V. SUMMARY 



A whole process of deriving the effective interaction in many-nucleon systems has been shown. What is presented 
in this chapter is a kind of modelling: the modelling of interacting many-nucleon systems. In particular we have 
illustrated the appearance of both nonlincarity and differential operators in this formalism. In the application of vari- 
ational principle, based on the functional analytic methods, we have presented a treatment of the Gateaux differential 
in some generalized situations. 

Phenomenological density-dependent force has been treated. The exchange part of the Coulomb force and the 
pairing force have fractional powers of the density, so that the density-dependent property is also true to these two 
forces. Such a fractional power dependence reasonably appears in the physics treating finite quantum systems. Indeed, 
as far as finite quantum systems are concerned, the existence of p term in energy implies the emergence of p 2 ^ 3 term 
in energy due to the surface effects. 

Two unknown features have been found. First, as is shown in the first term of the right hand side of Eq. (|16p . the 
coefficient of the kinetic e nerg y part of the Hamiltonian density is equal to h 2 /Am. This fact requires the modification 
of the standard choice [H, Il3l | of coefficients contained in the Hamiltonian density. Second, as is shown in Eq. (|33[) , 
we have pointed out that there is a missing term in the standard derivation of the effective interaction arising from 
j 2 in the Hamiltonian density. This fact requires the modification of the standard choice Q of coefficients contained 
in the effective interaction. The correct treatment of this term is necessary not only to hold the Galilean invariance, 
but also to clarify the time-odd contribution to the stationary and non-stationary states of many-nucleon system. 
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VI. APPENDIX -DIFFERENTIABILITY - 



Let z be a notation for complex variable: 

JL — 1 ( J> -iJt\ l-i/'l + .-l'l 

dz ~ 2 \dx 1 dy ) ' dz ~ 2 \ Ox ~ 1 dy J ' 

We consider a function f(z) satisfying f(z) = f(x, y) and z = x + iy. Here let us assume that f(x, y) is a differentiable 
function of the two real variables (x, y). According to the Taylor's theorem, 

f(z + u>) - f(z) = f{x + LV + T,)- fix, y) = + y%& V + ■ ■ • , 

where lu = £ + r\ and 



Therefore 



_s_ = _a_,_a_ d_ _ ^ t_d 

dx dz dz ' dy V dz dz > 



+ 



= s fia« + ft7) + 2 » i K -<»>) + ■•■ (37) 

= M£) W + M£)^ + .... 

dz dz 

In particular / is Frechct differentiable if df /dz = (i.e., / is holomorphic). However, even though / is not 
holomorphic, 

lim A ^ a - -dz-^ + -dT u 

is valid to real A, if / is differentiable. In this situation it is possible to consider Gateaux differential: 

nff y x dfiz) dfiz)_ 

oz az 

In particular if / is real-valued, the Gateaux differential is reduced to 

Df(z)(w) = 2Re (^w) , (38) 

where note that df(z)/dz and dfiz)/dz are complex-conjugate. In the case of Sec. ([Ill C 1|) . the right hand side of 
Eq. ([55)1 corresponds to 



-^SU) + \4>)\ x =o = -Re / dr 3 |Wr 2 V^V0, 



dX 



20 



but the limit (A — > 0) cannot exist if A is not real (on the other hand, A should be complex with respect to finding 
the optimal condition in the complex Banach spaces). Indeed, according to Eq. (|37p . the limit of 

ew+xtf-ew = dew , , MM at, ... 

A dip dip A™ 

due to A — > (A is not real) exists only when E is holomorphic. Note that £ shown in Sec. ()III C lj) is not holomorphic. 

In summary a mapping from a complex Banach space to a complex Banach space is Frechet differentiable, only 
if the function with complex variables is holomorphic. On the other hand, the real-valued functions are not holo- 
morphic except for the constant. It provides a reason why we have lengthy treatments of the Gateaux differential as 
demonstrated in Sec. IIII CI 



